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ABSTRACT 

In this paper, we have obtained existence, uniqueness and error bounds for deficient discrete quartic spline 
interpolation. 
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1. INTRODUCTION 

Let us consider a mesh on [0, 1] which is define by 

0=X 0 <X l < <*„=1 with X i —X i _ l =P i fori = 1,2,... .n. 

and h > 0, will be a real number, consider a real continuous function s(x,h) defined over [0, 1] which is such that 
its restriction S- on \x f _ x ,X J is a polynomial of degree 4 or less for i = 1, 2,....n. Then s(x, h) defines a discrete 
quartic spline if 

Dl l} s i {x i ,h)=Dl J] s i+l {x i ,h)j = 0,l (1.1) 
Where the difference operator D h are defined as 
Dl 0] f(x) = f(x) 

D u, m= f(*+h)-m) 

h 

D i2 } f{x + h)-2f(x) + f(x-h)) 

h 1 2h 

Let D(4-,l,A,h) is the class of deficient discrete quartic spline interpolation of deficiency one, where in 
D *(4,1, A,h) denotes the class of all discrete deficient quartic splines which satisfies the boundary condition 

s(x 0 ,h)=f(x 0 ,h) 

s(x n ,h)=f(x n ,h) (1.2) 
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Mangasarian and Schumaker [6, 7] introduced discrete quartic splines to find minimization problem. Existence, 
uniqueness and convergence properties of discrete cubic spline interpolation matching the given function at mesh point 
have been studied by Lyche [4, 5] which have been generalized by Dikshit and Power [1] (see also Dikshit and Rana [2]). 
It has been by Baneva, Kendall and Stefanov [3] that the local behaviour of the derivative of a cubic spline interpolator is 
some times used to smooth a histogram which has been estimated in [8], Rana [8] has obtained a precise estimate 
concerning the discrete cubic spline interpolating the given function at the mesh points (See also [9], [10]). In the present 
paper we obtain a similar precise estimate concerning the deficient quartic spline interpolant matching the given function at 
two intermediate points between successive mesh points and first difference of interior points of interval [0, 1]. 

2. EXISTENCE AND UNIQUENESS 

We introduced the following interpolating conditions for a given function/. 

s(flr / )=/(ar i ) 



D™s(r t )=D™f(y t ) 



(2.1) 
(2.2) 
(2.3) 



Where CC i =X i _ [ + ^P i = f j , for i = 1, 2,....n. 

we shall prove the following. 
Theorem 2.1 

Let f be a 1 -periodic, then for any h>0 then these exist a unique 1 -periodic deficient discrete quartic spline s in the 
class D*(4,l,A,h) which satisfies interpolatory condition (2.1) - (2.3). 

Proof 



Let P(z) be a quartic polynomial on [0, 1], then we can show that 

z) + Dl li p(\)q,(z) + P(0) 
V-V 





'p 




^) 


p 


■q l (z) + P 






v2j 



(2.4) 



Where q l (z) = 



1 3, 2 ^ L 9 ^ 
+— h 



v 6 2 



z+ 



2--h 2 
2 J 



2 . 3 

z +z 



29 



+ 24h 2 ) + (3-\Sh 2 )z 4 ] 



'21 ^ 
---h 2 
81 3 
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32 
81 



-z + 



f -224 16,, 

+ — h 

81 3 



2 . 3 

z +z 


O60 


64 

— h 


+ z 4 


f 


— +\ee 


J 


V 27 


3 j 






9 J 



f 2 1 ^ 



81 3 



J 



-1 2 2 11 3 2 4 

ZH — Z Z + — Z 

9 3 9 3 



'21 ^ 



81 3 j 



1 + 



V 



2 4 , , 
27 3 



z + 



J 



58 /? 2 ^ , f 26 16, 2 ^ 

H ft 



V 



81 + 3 



27 



z 3 + 



+ 4h 2 



81 3 



2 A 



1 j. 1 ^ 

h — ft 

162 6 



z + 



J 



162 6 



2 , 

z + 



7 8, 2 
54 3 



z 3 + 



1 



— 2h 2 
V9 j 



f 2__h^ 
81 3 



X-X ; 



Denoting t = -,0<?<1 , we can expressed (2.4) in the form of restriction s(x, h) of the deficient discrete 

Pi 

quartic spline s(x, h) on \x M , X M ] as follows :- 

S(x,h) = f {a, ) q x (x) + f (ft ) q 2 (x)+ P, D' 1} (ft ) (x) 

+ S i 0O# 4 (x) + S w (x)<? 5 (x) (2.5) 

Observing (2.5) it may easily be verify that S t (x, h)is a quartic on [3C-,JC i+1 ] for i=0, 1,...., n-1 satisfying 
(2. 1) - (2.3) and writing H(a,b) = Cl + bh , for real a, b we shall apply continuity of the first difference of s(x, h) at 
X: in (2.5) to see that 



H 



92 



81 



,2 



J 



2326 -32 



189 



i* 2 k 



+ 



\pAh 



-4 13 

v 27'6y 



-17 16 



54 3 



+ P.-AH 



< -2 4 
9"* 3 



v 



P-H 



26 16 
27' 3 
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+ P.-AH 



= F : 



1 1 



P; +H 



162 6, 

i= 1,2, n-1. 



Where F x =P i 



H 



-37 



59 



-7 8 
54 '3 



h 2 \s 



M 



(2.6) 



P 2 x +h 2 H 



~ 569 -48*0 

J 



fia^-P^H 
MJ-/(fl)|- 



-1 -3 

6 ' 2 



P 2 +H 



■29 



V 



V 6 
) 



.24 



f 128 


736 A 




I 81 


9 , 


V 



-224 128 V 7 



f 32 
H—,0 
81 



V 



27 3 



160 -64 

27 ' 3 



r 2 13) 
v9' 9 j 



V 9'9y 



Existence, uniqueness of s(x, h) depend on the existence of a unique solution of set of equation (2.6). It is easy to 
observe that in (2.6) absolute value of the coefficient of s, dominates over the sum of the absolute values of the coefficient 

of S M and S t _ { i.e. is positive. 



TXP,h) = 



H 



^80 O 
81 '6 



p f-i + H 



0511 16\ 2 TT f 35 7 
h +H 



126 3 



162 6 



P ; 2 +H 



45 8 

V 54'3y 



h 2 



Thus the coefficient matrix of equation (2.6) is diagonally dominant and hence invertible. 

Remark: In the case h— >0 theorem 2.1 gives the corresponding result for continuous quartic spline 
interpolation under condition (2.1) - (2.3). 

3. ERROR BOUNDS 

Now system of equation (2.8) may be written as 

A(h) M(h)=F 

Where A(h) is coefficient matrix and M{h) = S j {K) . However, as already shown in the proof of theorem 2.1. 

A(h) is invertible. Denoting the inverse of A(h) by A~'(h) we note that row max norm A~'(h) satisfies the following 
inequality :- 



\\A- l (h)\\<y(h) 



(3.1) 
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Where j(/l) = max{7^ (h)} ' . For convenience we assure in this section that 1 = Nh when N is positive integer. 

It is also assume that the mesh points { X- t } are such that X i G [0,1]^ for i = 0, 1 n. Where discrete interval [0,l] A is the 

set of points {0,h Nh} for a function/and two distinct points X x , X 2 in it domain the first difference is defined by 

Xj x 2 

For convenience, we write f for Dj}f and w(f , p) for modules of continuity off, the discrete norms of a 
function /over the interval [0, l] h is defined by 

II / 11= max I f(x) I (3.3) 

.re [0,1],, 

We shall obtain in the following the bound of error function e(x) = s(x, h) — f(x) over the discrete 
interval [0,1] h . 

Theorem 3.1 

Suppose s(x,h) is the discrete quartic spline interpolant of Theorem 2.1. Then 

\\e(x)\\<k(P,h)w(f,p) (3.4) 
\\e(x i )\\<y(h)K*(P,h)w(f,p) (3.5) 

We^ix^K^P^wif^) (3.6) 
Where K (p, h), K*(P,h) and K^P,^ are some positive functions of p and h. 

Proof 

Writing f(x i ) = f i equation 3.1 may be written as 

A(/l).( e (x,))=F,(/l)-A(/l)(/,.) = L,.(/) (Say) (3.7) 

When we replace S^h) by e i (x i ) = s(x i ,h)- f i (3.8) 

We need the following Lemma due to Lyche [4, 5], to estimate inequality (3.5). 
Lemma 3.1 

Let {a ; }™j and }" _^ be given sequence of non negative real numbers such that 

m n 

i=i j=i 
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Then for any real valued function/, defined on discrete interval [0, l] h , we have 

2>- k >**> • h I ~i b J k ' y 3i y h \f 



H 



<w(f ik \\l-Kh\) 



(3.9) 



Where X ■ ,y , e[0,l] A for relevant value of i, j and k. We can write the equation (3.9) is of the form of error 



h ' J h 
function as follows 



H 



(92 9 ) 



2326 -32 



+ 



r 



P \H 



■4 13 



V 



27 6 

f 2 ^ 



P^+H 



189 3 j 
r \l 16^ 



h 2 \e 



i— 1 



V 



54 3 



h 1 



+ P 



H 



V 9 3y 
' 1 1 



Pf-H 



26 16 
27' 3 



v 



162 6 



P 2 +// 



'-7 
v54'3y 



e M =R,(f) 



Where R.(f)=F-P? 
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(3.10) 



Writing equation (3.9) is the form of divided different and using Lemma 3.1 given by Lyche [5] 



|«,(/)|=h 3 ^k,A4]Jfl 
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(3.12) 
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b =P P 

u \ r \ r i-l 



H 



H 
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J 
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b 2 =P x P^H 



11 -1 



81 3 



h 1 3 3 

h = -gPiP i -l 



K=pliPi 



H 
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81 9 
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n 



81 9 



b 5=-j h2 pliPi 



and x 1q = = y h , x h =x,= x 2] = y 4] 



■^2 n i 



Ji 0 = «h . y 2o = ®i = a = *4 0 = * 5o 

y3 0 =7i-h = y 5o ,y 3i =y- t +h= 5i 

Now using the equation (3.10) and (3.9) in (3.8). 

\\e( Xi )\\<y(h)K*(P,h)h(f w ,p) 

This is the inequality (3.5) of Theorem 3.1. 

To obtain inequality (3.4) of Theorem 3.1. Writing equation (2.5) in the form of error function as follows: 
e(x)=e i _ 1 Q 4 (t) + e l Q 5 (t)+M i (f) 

Where M l (f ) = / (a, ) ft (f ) + / (# ) ft (0 



(3.13) 
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Denoting t = -,0<?<1 , we can expressed (2.4) in the form of restriction s(x, h) of the deficient discrete 

Pi 

quartic spline s(x, h) on \x M , X M ] as follows :- 

S(x,h) = f {a, ) q x (x) + f (ft ) q 2 (x)+ P, D' 1} (ft ) (x) 

+ S i 0O# 4 (x) + S w (x)<? 5 (x) (2.5) 

Observing (2.5) it may easily be verify that S t (x, h)is a quartic on [3C-,JC i+1 ] for i=0, 1,...., n-1 satisfying 
(2. 1) - (2.3) and writing H(a,b) = Cl + bh , for real a, b we shall apply continuity of the first difference of s(x, h) at 
X: in (2.5) to see that 
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+ PJ>™ftr t ) G 3 W (* ) + s, (x) Q { 4 li (t) + s i+1 Q? ] (t) 
A i P i e m (x) = e t _ x Qf ] (t) + e, (t) + U t (f) 

Where U i (/) = /(*.) (t) + /( ft ) Qf (f) + f{y { ) Q,f ] (0 

+ / w Q.t ) {t) + f i QV{t)-AP i f m {x) 
By using Lemma 3.1, we get 
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X l 0 = X i' X l l = Pi = ^l, 



(3.15) 
(3.16) 



(3.17) 
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y 2o =x-h,y 2j =x + h 

By using (3.5), (3.13) and (3.17) in (3.16), we get inequality (3.6) of Theorem 3.1. 
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